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ABSTRACT 

MHD turbulence consists of waves that propagate along magnetic fieldlines, in both 
directions. When two oppositely directed waves collide, they distort each other, without 
changing their respective energies. In weak MHD turbulence, a given wave suffers 
many collisions before cascading. "Imbalance" means that more energy is going in one 
direction than the other. In general, MHD turbulence is imbalanced. A number of 
complications arise for the imbalanced cascade that are unimportant for the balanced 
one. 

We solve weak MHD turbulence that is imbalanced. Of crucial importance is that 
the energies going in both directions are forced to equalize at the dissipation scale. We 
call this the "pinning" of the energy spectra. It affects the entire inertial range. 

Weak MHD turbulence is particularly interesting because perturbation theory is 
applicable. Hence it can be described with a simple kinetic equation. Galtier et al. 
(2000) derived this kinetic equation. We present a simpler, more physical derivation, 
based on the picture of colliding wavepackets. In the process, we clarify the role of the 
zero-frequency mode. We also explain why Goldreich & Sridhar claimed that perturba- 
tion theory is inapplicable, and why this claim is wrong. (Our "weak" is equivalent to 
Goldreich & Sridhar's "intermediate.") 

We perform numerical simulations of the kinetic equation to verify our claims. We 
construct simplified model equations that illustrate the main effects. Finally, we show 
that a large magnetic Prandtl number does not have a significant effect, and that hy- 
perviscosity leads to a pronounced bottleneck effect. 

Subject headings: MHD — turbulence 



1. Introduction 

MHD turbulence is ubiquitous in astrophysics. For example, it is present in the sun, the solar 
wind, the interstellar medium, molecular clouds, accretion disks, and galaxy clusters. Theoret- 
ical understanding of incompressible MHD turbulence has grown explosively in the last decade. 
Nonetheless, it remains underdeveloped. 

Iroshnikov (1963) and Kraichnan (1965) developed a theory for MHD turbulence. They realized 
that the magnetic field at the largest lengthscale in a cascade directly affects all of the smaller 
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lengthscales. Small-scale fluctuations can be treated as small-amplitude waves in the presence of 
a large mean magnetic field. By contrast, the large-scale velocity is unimportant for small-scale 
dynamics; it can be eliminated by a change of variables, since the equations of MHD are invariant 
under Galilean transformations. 1 

Despite their realization of the importance of the mean magnetic field, Iroshnikov and Kraich- 
nan assumed that small-scale fluctuations are isotropic. Numerical simulations later showed that 
this assumption is wrong. Even with isotropic excitation at large scales, fluctuations on smaller 
scales are elongated along the mean magnetic field (e.g., Montgomery & Turner 1981, Shebalin, 
Matthaeus, & Montgomery 1983). 

In retrospect, it is not very surprising to find elongated fluctuations. Arbitrary disturbances in 
incompressible MHD can be decomposed into shear-Alfven and pseudo-Alfven waves. Each wave 
travels either up or down the mean field at the Alfven speed, va, which is the magnitude of the 
mean field in velocity units. Consider stirring a magneto-fluid with a spoon that is moving at speed 
v <C va, for a time comparable to the spoon's width divided by v. 2 Alfven waves are radiated away 
from the spoon, parallel to the mean field with speed ±va- After the disturbance is finished, there 
are two wavepackets travelling away from each other. Each wavepacket is elongated along the mean 
field, with parallel-to-transverse aspect ratio ~ va/v S> 1. 

The characteristics of MHD turbulence depend critically on the amount of elongation. When 
parallel-to-transverse aspect ratios are smaller than va/v, waves collide many times before cas- 
cading. Hence the turbulence is "weak," and perturbation theory can be used to derive a kinetic 
equation and a spectrum (Sridhar & Goldreich 1994; see Zakharov, L'vov, & Falkovich 1992 for a 
general review of weak turbulence). Goldreich & Sridhar (1997) and Ng & Bhattacharjee (1997) 
deduced the spectrum of the balanced weak cascade from scaling arguments. 3 However, similar 



x One outgrowth of Iroshnikov and Kraichnan's theory that is particularly relevant to the present paper is Grappin, 
Pouquet, and Leorat (1983). We discuss it in §4. 

2 In a turbulent cascade, one would expect v <C va on small scales, since v decreases towards small scales, whereas 
va is unchanged. 

3 We relegate some of the history to a footnote because it can be confusing. Sridhar & Goldreich (1994) developed 
the first theory of MHD turbulence that accounted for the anisotropy of fluctuations. They claimed that "three-wave" 
processes vanish in weak MHD turbulence, and "four-wave" processes must be considered (i.e., perturbation theory 
is trivial to first order, so second order terms are important). As a result, they used four- wave couplings to derive a 
kinetic equation and a spectrum for weak MHD turbulence. Montgomery & Matthaeus (1995) claimed, and Ng & 
Bhattacharjee (1996) showed, that Sridhar & Goldreich (1994) are wrong, and three-wave processes do not vanish. 
Goldreich & Sridhar (1997) explained the contradiction: Sridhar & Goldreich (1994) had unknowingly assumed that 
fieldline wander is limited, i.e., that the separation between any two fieldlines is nearly constant along their entire 
length; in this case, three-wave couplings are negligible and the kinetic equation based on four-wave couplings is 
correct. In the more realistic case that fieldlines do wander, three-wave processes are important. Goldreich & Sridhar 
(1997) went on to argue that, in the latter case, perturbation theory is inappropriate, and couplings of all order are of 
comparable magnitude; so they called this "intermediate" turbulence. Galtier et al. (2000) argued that perturbation 
theory is appropriate, even when three-wave processes are important. In the Appendix of the present paper, we use 
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scaling arguments are inadequate for the imbalanced cascade (see §3.2 of the present paper). Galtier 
et al. (2000) derived the kinetic equation for the weak imbalanced cascade. Their balanced spec- 
trum agrees with that of Goldreich & Sridhar (1997) and Ng &: Bhattacharjee (1997). They also 
presented a partial solution for the general imbalanced case. In §3.2, we explain why their solution 
is incomplete; in §4 we give the complete solution. 

Even if aspect ratios are smaller than va/v on large scales, at a small enough scale they become 
comparable to va/v. Below this scale perturbation theory breaks down, and weak turbulence 
becomes "strong." Goldreich h Sridhar (1995) worked out the scalings for the balanced strong 
cascade. They argued that aspect ratios are comparable to va/v at all scales in the strong regime. 
Strong turbulence is difficult, largely because it is non-perturbative. Although strong and weak 
turbulence differ in a number of ways, they also share many similarities. One of our motivations for 
studying weak turbulence is to gain insight into strong turbulence. In particular, turbulence in the 
solar wind is observed to be imbalanced; it cannot be understood without a theory for imbalanced 
strong MHD turbulence. Yet this theory is unknown. In a future paper, we will work it out by 
extending the results of the present paper. 



2. Basic Equations 

Ideal incompressible MHD 4 is described by the following equations of motion: 

dtV + v-Vv = -VP+B VB, (1) 
d t B + v-VB = B Vv , (2) 
V • v = V • B = . (3) 

The density is set to unity; the fluid velocity is v; the magnetic field in velocity units is B = (magnetic 
field)/ (4-7T) 1 / 2 ; the total pressure is P = p + B 2 /2, which is the sum of the thermal and magnetic 
pressures. Viscous and resistive terms are neglected in the above equations; they are important on 
small scales, and will be included where required. 

We decompose the magnetic field into its mean, vaz, where va is the Alfven speed and £ is a 
unit vector, and into its fluctuating part b = B — vaz. With this decomposition, the equations of 
motion may be written in terms of the Elsasser variables, = v — b and = v + b, as follows: 

d t w^ + v a9 z w^ = -w^-Vw^ - VP , (4) 
dtw^ — va9 z w^ = —w^ -VwJ- — VP , (5) 
V-io T = V-tu* = . (6) 



Goldreich & Sridhar's picture of wavepackets following wandering fieldlines to clarify the controversy, and to explain 
why perturbation theory works. Because it does work, we call the cascade "weak" instead of "intermediate." 

4 In this paper, we consider only incompressible MHD turbulence; compressibility does not alter the dynamics very 
much (Lithwick & Goldreich 2001). 
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Note that P is not an independent degree of freedom. Taking the divergence of either equation (4) 
or (5) yields 

P = -V" 2 (Vi« t :Vw±) , (7) 

where V -2 is the inverse Laplacian. When = 0, propagates undistorted upwards along the 
mean magnetic field with speed va- Similarly, when = 0, propagates downwards at va- 
Nonlinear interactions occur only between oppositely directed wavepackets. It is these interactions 
that are responsible for turbulence. 

There are three conserved quantities in incompressible MHD. Two of these are immediately 
apparent from equations (4)-(6): the energies of the upgoing and of the downgoing waves, i.e., 
(w 1 *) 2 and (w^) 2 . Technically, these are twice the energy per unit mass. We refer to them as 
simply energies throughout the paper. These energies are directly related to the total (kinetic plus 
magnetic) energy oc (w 1 ') 2 + (w^) 2 and to the cross-helicity oc (w^) 2 — (w^) 2 . The focus of this 
paper is turbulence where the energies in the up and down waves differ, or, equivalently, where the 
cross-helicity is non-zero. The third conserved quantity is magnetic helicity; however, we consider 
only non-helical turbulence in this paper, so helicity does not play a role. 

In MHD turbulence, on lengthscales much smaller than the outer scale, there is effectively 
a strong mean magnetic field that is due to fluctuations on the largest lengthscales. Gradients 
transverse to this mean field are much larger than gradients along it (e.g., Shebalin, Matthaeus, 
& Montgomery 1983, Goldreich & Sridhar 1995, 1997, Ng & Bhattacharjee 1996). This allows 
the MHD equations to be slightly simplified. Denoting transverse components with the symbol 
_L = (x,y), the transverse components of equation (4) are 

d t w]_ + v A d z w]_ « -w\_- V ± w]_ - V ± P , (8) 

assuming that w z d z w^ is much smaller than w^_- Vj_tu^ , and equation (6) is 

V±-w]_ wO. (9) 

We assume that the parallel components of and are either comparable to, or less than, 
their respective perpendicular components. We will see below that this is typically the case in the 
inertial range of a turbulent cascade. Similarly, 

d t w l ± - v A d z w l ± « -w]_ -V_li<;^ - V_lP , (10) 



Vi^RiO. (11) 

If we change « to =, equations (8)-(ll) form a closed set. They are called the equations of reduced 
MHD. They apply also in compressible MHD whenever transverse gradients are larger than parallel 
ones (e.g., Biskamp 1993). The main goal of this paper is to solve these equations. Although the 
complete equations are not much more complicated, it simplifies our discussions to neglect parallel 
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gradients relative to perpendicular ones at the outset. There are two conserved energies in reduced 
MHD: (u>|) 2 and {w l ± ) 2 . 

The parallel components of equations (4) and (5) are 

d t w\ + v A d z w\ « -w l ± -V ± wl , (12) 
dtw^ — va9 z w\, — -Viw] , (13) 

T i ili 

after neglecting parallel gradients relative to transverse ones, and after assuming that \w z \ and \w z \ 
are not much smaller than \w\_ | and Clearly, (wl) 2 and (wl) 2 are conserved quantities. The 

transverse equations describing reduced MHD are unaffected by these parallel equations because 
the former are independent of wl and w\. Nonetheless, the parallel equations have observable 
consequences. 

It is conventional to decompose the normal modes of linearized incompressible MHD, and 
w^, into shear-Alfv'en and pseudo-Alfven waves. These correspond to the Alfven and slow waves 
of compressible MHD. When perpendicular gradients are much larger than parallel ones, w\_ and 
w\_ are nearly equivalent to shear- Alfven waves; wl and wl are nearly equivalent to pseudo-Alfven 
waves. 

Also observationally relevant is the evolution of a passive scalar s, which satisfies dts + v-'Vs = 
0. In terms of Elsasser variables, and after neglecting parallel gradients, the passive scalar satisfies 

d t s « -(l/2){w]_ + w l ± )-V ± s . (14) 

To avoid a proliferation of subscripts, in the remainder of this paper we drop the _L from w\_ 
and w\_. To denote the parallel components, we use wl and wl. 

3. Weak MHD Turbulence: Heuristic Discussion 

One of the virtues of weak MHD turbulence is that it can be analyzed in a mathematically 
rigorous way with perturbation theory; this yields a kinetic equation. Nevertheless, we begin with 
a qualitative description, which captures most of the features of the turbulent cascade. 

3.1. Scaling Relation 

MHD turbulence can be understood from the dynamics of and (eqs. [8-11] for reduced 
MHD, dropping _L subscripts). To linear order, is a wave that propagates up the mean magnetic 
fieldlines at the Alfven speed, va', propagates down at va- Each wave perturbs the mean 
magnetic fieldlines. Nonlinear terms describe the interaction between oppositely directed waves: 
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each wave nearly follows the fieldlines perturbed by its collision partner. 5 

Consider an upgoing wavepacket that encounters a train of downgoing wavepackets. As the 
upgoing wave travels up the length of the downgoing train, it is gradually distorted. It tries to follow 
the perturbed fieldlines in the downgoing train, but these fieldlines "wander," i.e., the transverse 
separation between any two fieldlines changes. When the up-wave has travelled through a suffi- 
ciently large number of downgoing wavepackets that the amount of fieldline wander is comparable 
to the up-wave's transverse size, then the up-wave cascades. 

To be quantitative, let each downgoing wave in the train have a typical amplitude w[, a 
transverse size A, and a parallel size A, where "transverse" and "parallel" refer to the orientation 
relative to the mean magnetic field. The most important collisions are between wavepackets of 
comparable transverse size (see §4.3.1). So let the upgoing wave have transverse size A as well. 

Since each downgoing wavepacket has a typical perturbed magnetic field of magnitude ~ 
(neglecting the factor of 1/2), it bends the fieldlines by the angle w^/va; the transverse displacement 
of a fieldline through this wavepacket is (w\ /va)A; and the wander of two typical fieldlines through 
the wavepacket, if they are initially separated by A, is also (w^/va)A. 

In weak turbulence, the wander through a single wavepacket is smaller than the wavepacket's 
transverse size, 

I T 
w w 

-^A < A and -^A < A . (15) 

va v A 

When these inequalities are not satisfied, strong turbulence is applicable; see §4.3.2. Thus in weak 
turbulence an upgoing wavepacket must travel through many downgoing ones before cascading. 
After N downgoing wavepackets, fieldlines have wandered a distance ~ N l / 2 (w ^ x /va)A, assuming 
that wavepackets are statistically independent. The upgoing wavepacket is fully distorted — and 
hence cascaded — when the fieldlines it is following wander a distance A, i.e., when N ~ (Xva/ Aw{) 2 . 
Since each downgoing wavepacket is crossed in the time A/va, the cascade time of the upgoing 
wavepacket is 



^~(^) 2 M4)x- < 16 > 

y Awy v A ^wy A 

In this time, the upgoing wavepacket travels a distance v Atlas, which is much larger than its 
own length, A. 6 The head and the tail of the upgoing wavepacket are both distorted by the same 



J The equation for a scalar quantity / that travels upwards at speed va, while following the magnetic fieldlines 
of the down-going is (dt + vaQz + w^-Vi)/ = 0. Equation (8) for the vector differs from this because 
of the pressure term, which is required to keep incompressible, while conserving the energy (w^) 2 . Thus 
does not exactly follow the fieldlines of io J \ Nonetheless, this deviation does not greatly affect the behaviour of 
the turbulence. Dissipation is a second effect that prevents the following of fieldlines. In the present discussion, we 
consider lengthscales that are sufficiently large that dissipation can be neglected. 

6 We assume throughout this paper that the upgoing waves' parallel lengthscale is the same as that of the downgoing 
waves, A; the extension to the case when they differ is trivial, as long as the inequalities (15) are both satisfied, with 
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downgoing wavepackets; so both head and tail undergo nearly the same distortion as they cascade. 7 
Consequently, as the upgoing wavepacket cascades to smaller transverse lengthscales, it does not 
cascade to smaller parallel ones: 

A = scale independent . (17) 

A proof of this follows from the 3-wave resonance relations (Shebalin, Matthaeus, Sz Montgomery 
1983; see also §A.5 of the present paper). 

We calculate the steady state energy spectra by using Kolmogorov's picture of energy flowing 
from large to small lengthscales. The energy in up- waves flows from lengthscales larger than A to 
those smaller than A at the rate 



h2 



e T^ K) 

' cas 



T 1 

w x w x 
A 



2 A 



(18) 

VA 



We call this simply the "flux." We define more precisely below (eq. [51]). In steady state, the 
flux must be independent of A, so 

w{w l x oc A . (19) 



3.2. Insufficiency of Scaling Arguments for the Imbalanced Cascade 

A balanced cascade has w\ = = w\. Its solution in steady state is simple: w\ oc 
and e T = e l ~ w{k/\ 2 v A (Goldreich & Sridhar 1997, Ng & Bhattacharjee 1997). However, if the 
cascade is imbalanced, a number of complications arise. 

By the symmetry between up- and down-going waves, the down-going flux is given by the 
analogue of equation (18): 



£-1 ~ 



T I 
w x w x 



A 



2 A 



(20) 

VA 



Because and both depend on the same combination of w\ and w x — namely their product — the 
steady state solution is non-trivial. Had this degeneracy not occurred, e.g., had we found 

T (w\) 2+ ^w{) 2 -^ A , , (^!) 2 " 7 M) 2+7 A , , 

£ T „ v A; V X! and e i „ K X) V a; 21 

\ z v A y V A 

where 7 7^ 0, then the solution would have been simple: w\ oc w oc A 1 / 2 and (eVe 1 ) ~ {w\/w[) 2 \ 
which follow from the constancy of and with A. 



the appropriate A's. 

7 The head of the upgoing wavepacket slightly distorts each downgoing wavepacket; so the downgoing wavepacket 
seen by the tail is slightly distorted relative to that seen by the head. Nonetheless, this backreaction is a higher-order 
correction that can be ignored in weak turbulence; see §A.l. 



-8- 



But in weak MHD turbulence 7 = (eqs. [18] and [20]); constancy of with A is forced by 
the constancy of e^, and does not yield new information. One implication is that scaling arguments 
are insufficient to determine the flux ratio /eK Physically, any flux ratio should be possible. 
But without the dimensionless coefficients of equations (18) and (20), e^/e^- cannot be determined. 
The coefficients cannot be obtained from scaling arguments; they depend on the spectral slopes of 
w\ and w x (which are related through w\w x oc A, eq. [19]). Galtier et al. (2000) calculated the 
coefficients using kinetic equations. We explain how in §A.4. Therefore these authors were able to 
relate the flux ratio to the spectral slopes. 

The arguments presented thus far are still insufficiently constraining. Equations (18) and (20) 
constrain only the product w^Wy There are seemingly an infinite number of solutions with given 
values of and e^, since w]^ can be multiplied by any constant as long as w x is divided by this 
same constant. Furthermore, we expect on physical grounds that if the values of w x and w x at a 
given lengthscale are fixed (instead of the values of and e^), the cascade should be completely 
constrained; however, in this case the constancy of equations (18) and (20) leaves the A-dependence 
of w\/w x completely undetermined — even given the coefficients derived by Galtier et al. (2000). 

T I 

Do wl and w\ cross? Are they cut off by dissipation at the same scale? All of these problems for 
the imbalanced cascade can be resolved once the dynamics at the dissipation scale is understood. 



3.3. Dynamics at the Dissipation Scale: Pinned Spectra 

The main result of the present paper is that the energies of the up- and down-going waves 
are forced to equalize — they are "pinned" — at the dissipation scale. This completely constrains the 
cascade. It is unusual that the dynamics at the dissipation scale has such an important influence. 
In this subsection we explain why pinning occurs. In §4, we give the resulting solution of the steady 
state cascade. 

From equation (16), the cascade time of the up-going waves is inversely proportional to the 
energy of the down-going ones: i<t a s = (^/w x ) 2 (va/A-), and similarly for the downgoing waves. We 
consider how the spectra evolve if initially, on lengthscales comparable to the dissipation scale, 
waves going in one direction are more energetic than the oppositely directed ones. To facilitate the 
discussion, we refer to Figure 1 on page 19, which presents the results from a numerical simulation 
that we discuss in detail in §6.1. In the middle panels of Figures la-d, we plot e} (k) ~ (Awlj 2 and 
e^(k) ~ (Au^) 2 as functions of wavenumber k = 1/A. The initial condition is shown in Figure la. 

t I T I 

Initially, w x > wt, so teas > ^cas- We consider a lengthscale-dependent dissipation time, tdiss> that 
is the same for both up- and down-going waves. 8 On large lengthscales, the dissipation timescale is 
much longer than the cascade times; t^iss decreases faster with increasing k than both tc as and teas- 
The effects of dissipation are felt on lengthscales where tdiss is comparable to — or smaller than — 



For example, in the simulation presented in Figure 1, i diss ~ A 2 /^, where v is both the viscosity and the resistivity. 
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either t cas or teas- Since teas > *cas in the vicinity of the dissipation scale, the largest lengthscale at 
which dissipation effects are felt is where teas ~ Miss- I n Figure la, this is at k ~ 4,000. We now 
let the spectra evolve; see Figs. lb-c. We hold the energies at k ~ 1 fixed; this does not affect the 
short-time behaviour shown in Figs. lb-c. Since w x feels the dissipative effects at k ~ 4, 000, its 
spectrum is exponentially cut off at smaller scales. This implies that the cascade time of the down 
waves , teas , increases exponentially towards smaller scales. As a result, down- wave energy that is 
being cascaded from large to small scales cannot be cascaded fast enough at k > 4, 000. Therefore 
the down-waves' energy flux is backed up, and the w[ spectrum rises. Furthermore, as w\ rises, 
teas falls, so the cascade time of the up- waves on small scales decreases, and the up-wave spectrum 
falls. The final result is that the two spectra are pinned at the dissipation scale. This pinning 
occurs very quickly: on the dissipation timescale. 



4. Steady State Energy Spectra 

T I 

In steady state, the energy spectra w x and w x are power laws that (i) are pinned at the 
dissipation scale, and (ii) satisfy w x w x oc A (eq. [19]). These two conditions completely characterize 

T I 

the steady state spectra, as long as w x and w x are specified at the outer scale. It is a remarkable 
feature of weak MHD turbulence that the dynamics at the dissipation scale dramatically affects 
the entire cascade. Normally, the energy in turbulent cascades is viewed as flowing unimpeded 
from larger scales to smaller scales. Energy is injected at the outer scale and swallowed up at the 
dissipation scale. But in weak MHD turbulence, if initially the spectra are not pinned, then the 
spectrum with lower energy becomes backed up, and its energy increases until pinning occurs. 

A similar effect was found by Grappin, Pouquet, and Leorat (1983). 9 These authors modelled 
imbalanced MHD turbulence with an EDQNM closure approximation. However, they assumed that 
the turbulence is isotropic, which is incorrect. Nonetheless, they found in their model that scaling 
arguments constrained only the sum of the slopes of the up- and down-going waves' spectra, and 
that the two spectra are constrained to be equal at the dissipation scale. Thus they discovered 
pinning, even though they analyzed an invalid model of MHD turbulence. 



In the remainder of this section, we derive the scalings in steady state. We denote the dissi- 

•r x and w x 



pation scale by Adi SS , and the value of wl and w x at Adi ss by Wdiss- We can then express the energy 



spectra as follows 

t / A \(l+a)/2 

w\ = w di J- (22) 

v -Miss 7 

. / A \(l-a)/2 , ^ 

w{ = w di J- . (23) 

v -Miss 7 

These spectra are valid in the "inertial range," i.e., on lengthscales larger than the dissipation scale, 



9 We thank Bill Matthaeus for pointing out this reference to us. 
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and smaller than the outer scale. There are three parameters that must be calculated to constrain 
the spectra: u>diss> Adi ss > and a. 

For definiteness, we assume that dissipation is caused by a diffusive process, described by a 
viscous term of the form u'V 2 v in equation (1) and a resistive term of the form vV 2 B in equation 
(2). This implies that the magnetic Prandtl number is equal to one. 10 The dissipative timescale is 

WA) ^ ^ • (24) 

At the dissipation scale, the cascade time is equal to t^ iss , i.e., £ C as(Adiss) = £diss(Adiss)> which implies 
that 




after using equations (16) and (24); so 

Wdiss ^ ("X") 1 ■ ( 26 ) 

Thus when the dissipation is caused by a diffusive process, Wdiss is independent of the outer-scale 
energies and the inertial-range fluxes. This is not true for Adi ss or a. 

To calculate Adiss and a, we consider two alternative scenarios: specified energies at the outer 
scale, and specified fluxes. 



4.1. Fixed Energies at the Outer Scale 

Suppose that the energies are specified at the outer scale, A ou tj where w\ and w\ are denoted 

^ ^A out aIld SillCe W l W \/ X = ^diss 2 /Adiss, 

. Aout vv A 

Adiss - —+ t 7 , (27) 

^out ^out 

after using equation (26). 

Dividing equation (22) by equation (23) yields w \ out / w \ out = (A ut/Adiss) a , so 

lnfu;! /w[ 1 In \w\ lw[ 1 

L Aout' A Qu t J ^ L A ou t / A ou t J (28) 

In [Aout/ Adiss] ln[w{ ou w l Xou A/(isv A )} 

With w\ / w \ out fixed, the cascade is balanced (a — > 0) in the limit that the inertial range is 
infinitely large (A ou t/ Adiss ~~ * 00 )• 



10 It is straightforward to consider other forms for the dissipation. We consider magnetic Prandtl numbers greater 
than unity in §8, and hyperviscosity in §9. 
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Inserting equations (26), (27), and (28) into the spectra (eqs. [22] and [23]) gives the solution 
to the steady state imbalanced weak cascade, assuming that w\ and ^j Vout are specified. 

Although we have solved for the spectra, recall from §3.2 that heuristic arguments are insuffi- 
cient for calculating the ratio of the fluxes that are carried by these spectra, e^/e^. In the following 
subsection, we show how ^ je 1 is related to the spectra. This relation is particularly important for 
solving the inverse problem: given the fluxes and e^-, what are the spectra of w\ and w^? 



4.2. Fixed Fluxes 



Physically, we expect that when the fluxes are fixed, the cascade should be completely con- 
strained. In this subsection, we solve for the spectra given and eK To accomplish this, the 
dimensionless coefficients of equations (18) and (20) are required. For a given power-law solution, 



w 



I oc \( 1+c *y 2 and w[ oc A^ 1 a ^ 2 , these coefficients depend on a: 



J - 



/(«) 
/(-« 



T I 
w x w x 



T I 
A 



2 A 

VA 
2 A 
VA 



(29) 
(30) 



where f(a) is a dimensionless function of a that must be calculated. By the symmetry between up 
and down waves, and are both proportional to the same function /, evaluated at ±a. Since 
heuristic arguments are insufficient to calculate the function /, it is fortunate that weak turbulence 
can be analyzed with perturbation theory. Galtier et al. (2000) computed /. We compute it in 
equation (A29) in the Appendix. 



The ratio of the fluxes is related to a by 



€' 

7i 



/(«) 
/(-«) 



(31) 



The limit \a\ <C 1 is particularly interesting. For given outer-scale energies, if the inertial range is 
very large then the steady state cascade is nearly balanced (see the discussion below eq. [28]), and 
\a\ <C 1. In this limit, we show in the Appendix that f(a) ~ /(0) • (1 + 0.5a) (see eq. [A30]). Thus 



e 1 

7i 



1 ~ a 



\a\ < 1 . 

To linear order in a, the product of the fluxes is independent of a: 



~ [f(0)f 



T I 



A 



^diss 



- -^diss 



(-Y 



\a\ < 1 



(32) 



(33) 



In the Appendix, we show that /(0) = 1.87 (eq. A31]). Although for most purposes the precise 
value of /(0) is unimportant, we shall need it when discussing our numerical simulations. 
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To summarize, in the limit of small \a\, if and are specified, then the spectra of and 
are given by equations (22) and (23), with w^iss, ct, and Adi ss given by equations (26), (32), 
and (33). Note that, to first order in a, the only relation that depends on the kinetic equation is 
equation (32). 11 

If, instead of specifying and , we specify the outer scale energies — in which case the spectra 
are given in §4.1 — then equations (32) and (33) give the resulting fluxes. 

4.3. Three Peripheral Issues 

This subsection may be skipped on a first reading, as it does not impact our main line of 
argument. 



4-3.1. Locality 

In §3.1, it is assumed that the dominant interactions are those between wavepackets that have 
comparable transverse lengthscales, i.e., interactions are "local" in lengthscale. In this section, we 
justify this assumption. 

We focus on the cascading of an upgoing wavepacket by downgoing ones. Let the upgoing 
wavepacket have transverse size A, and let the downgoing wavepackets each have transverse size 
I, parallel size A, and amplitude wj. The upgoing wavepacket cascades when the fieldlines it is 
following wander a distance comparable to its transverse size, A. 

We consider first the case that I < A. Two fieldlines that are separated by A at the head of the 
downgoing wavepackets wander independently of each other. Their transverse separation after N 
downgoing wavepackets increases by 

N 1 / 2 {w\/v A )K , I < A . (34) 

Conversely, if I > A, then the magnetic field can be expanded to linear order in A, so the 
magnetic field at two points separated by A differs by ~ w\\/l. Consequently, the fieldlines separate 
by 

N 1 / 2 (wj;/v A )(X/l)A , I > A , (35) 

so long as this separation is smaller than A. 

For interactions to be local, i.e., for the amount of fieldline wander seen by an up- wave of 
transverse size A to be maximized by those down-waves that have I ~ A, the following two conditions 



11 Aside from the uninteresting dependence of equation (33) on /(0). 
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must hold: (i) w\ is an increasing function of I (eq. [34]); and (ii) wj/l is a decreasing function of 
I (eq. [35]). So the cascade is local if 

dlnwl , . 

0< -5Ea <1 - (36> 

The same condition clearly holds for w\. In terms of the steady-state scalings, w\ oc X^+^l 2 and 
w[ oc \^~ a )/ 2 ^ the condition 

-1< a < 1 (37) 

is required for the cascade to be local; otherwise, nonlocal effects are important. Galtier et al. 
(2000) derived the inequalities in (37) from their kinetic equation. 

There is a second reason why the condition din wl /din A < 1 is necessary for our heuristic 
arguments to be valid. Consider a single Fourier mode with wavelength I and amplitude wl ; then 
the typical difference in w 1 * between two points separated by A < / is ~ wjx/l, expanding to linear 
order in A. So, if the condition din w\/ din A < 1 is violated, then the contribution to w\ (i.e., to the 
typical difference in between two points separated by A) is dominated by those Fourier modes 
that have wavelengths I > A. This is a different kind of non-locality than considered previously: 
the upgoing energy that crosses lengthscale A comes directly from much larger scales (with I S> A). 



4-3.2. Transition to Strong Turbulence 

Weak turbulence is applicable when w^A/va <C A and w^A/va <C A (eq. [15]). Since A 
decreases faster than both w\ and wl (eq. [36]), even if these inequalities are satisfied at large 
lengthscales, they are violated at small ones. Thus weak turbulence has a limited inertial range. 

Strong turbulence, which is applicable when the above inequalities are violated, is of greater 
relevance than weak turbulence for describing astrophysical sites such as the solar wind. In strong 
turbulence the change in fieldline separation within a single wavepacket is not smaller than the 
transverse size of the wavepacket. This has two implications. First, equation (16) for the cascade 
time is no longer valid; and second, the parallel size of wavepackets A decreases towards smaller 
scales because the head and tail of a wavepacket are independently cascaded. The balanced strong 
cascade is worked out by Goldreich & Sridhar (1995); we discuss the imbalanced strong cascade in 
a future paper. Strong turbulence is more difficult to analyze than weak turbulence because it does 
not submit to perturbation theory. Nonetheless, a number of the features of weak turbulence that 
we develop in the present paper are applicable to strong turbulence. This is one of our motivations 
for studying the weak cascade. 

Since this paper is concerned with weak turbulence, we choose the dissipation scale to be 
sufficiently large that the entire cascade is weak, i.e., we require that 

wl A 

A f s < 1 , (38) 
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and similarly for . Using equations (26) and (27), we can re-write this condition as follows: 



v\ A \Av A 

There is also a lower limit on the product of the outer scale energies, set by the requirement that 
the dissipation scale be smaller than the outer scale, which implies that 

T I 

w x ou ^x out > ^ (40) 



(see eq. [27]). Typically, one expects that A ~ A out , so as long as v <C Ava, the above two 
inequalities can be satisfied with appropriately chosen ^lout^Aout- 



4-3.3. The Steady State Spectra of w z and w z and of a Passive Scalar 

t i 

In §2 we showed that the parallel components w z and w z do not affect the evolution of the 
perpendicular components, and w^. 12 By contrast, the perpendicular components control the 
evolution of the parallel ones. 

The steady state spectra of the parallel components can be derived as follows. By analogy 
between the equation for w\ (eq. [13]) and that for (eq. [8]), the cascade time of w\ is similar 
to that of , i.e., it is given by i<t as (see eq. [16]). We denote the typical amplitude of the parallel 
component of the up- wave on lengthscale A by w\ x . The cascade rate of the energy of the parallel 
component is (w^ A ) 2 /tJas; in steady state it must be independent of A. Comparing this to the 
cascade rate of the perpendicular component, (w|) 2 /tj as , which must also be independent of A, we 
deduce 

w l,x « w x ■ ( 41 ) 

Note that the evolution equation for w\ (eq. [13]) is linear. So the overall amplitude of the w\ x 
spectrum is arbitrary; more accurately, it is set by the value of w\ A at the outer scale. Of course, 
a similar relation holds for the down- waves: 

w z,x « w x ■ ( 42 ) 

Finally, we consider the spectrum of a passive scalar, s, whose evolution is described by equation 
(14). From this equation, it is apparent that the cascade time of the passive scalar is given by either 
tj as or tc as , whichever is shorter. So, by the same reasoning that we used in deriving equation (41), 
we deduce 

sx oc mm(w{,w l x ) , (43) 



Recall that we drop the _L label. 
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where s\ is the typical value of s on lengthscale A. Because of the pinning of the spectra, the lesser 
of the two spectra w\ and is also the flatter. So s\ is proportional to the flatter of w\ and w^. 



5. Kinetic Equations in Weak Turbulence 

Weak turbulence can be analyzed with perturbation theory. As a consequence, the evolution 
of the energy spectra of the up- and down-waves is described by a closed set of two equations; in 
other words, the two-point correlation functions evolve independently of all higher-order correlation 
functions. This is a great simplification. 

Evolution equations for the energy spectra — the "kinetic equations" — were obtained in Galtier 
et al. (2000, 2001). We present an alternate, more physical, derivation in the Appendix. Such 
a derivation is useful because it clears up a number of erroneous claims that have been made in 
the literature — in particular, the claim of Goldreich & Sridhar (1997) that perturbation theory is 
inapplicable. 

In the following, we summarize the result derived in the Appendix. The kinetic equations are 
given in Fourier-space. We Fourier transform (x,y,z,t) and w^(x, y, z, t) in x and y (but not 
z), and denote the transforms by w^(z, t) and w^(z, t), where k is purely transverse (k z = 0). We 
define the energy spectra e' and e 1 such that 



*).«£,(*,*)) = e\k,t)5(k + k>) (44) 
w l k (z,t)-w l kl (z,t)) = e l (k,t)5(k + k') , (45) 



where 6(k) is a two-dimensional Dirac delta-function, and angled brackets denote both an ensemble 
average and an average over z. (We assume that the turbulence is homogeneous in z.) Both e} and 
are real. The turbulence is isotropic in the transverse plane, so e 1 * and e 1 are functions of the 
magnitude of k. 



d_ 

at 



From equation (A24) in the Appendix, the kinetic equation for the up-waves is 

e T (fc,t) = —k 2 [ dk 2 k^\e^{k 2 ,t)-e ] (k,t)} f dflsin 2 flcos 2 -uk 2 e\k,t) (46) 

V A Jo L -I Jo fci 



k 

where 

fci = (k 2 + kl- 2kk 2 cos 9) 1/2 . (47) 

We have included a term describing diffusive dissipation on small lengthscales, —vk 2 e\ where v is 
the viscosity, which is assumed to be equal to the resistivity. 

Because of the symmetry between up- and down-waves, the equation for e 1 is the same as that 
for e\ but with f and [ everywhere switched. The steady-state relation between the energy and 
the flux that we use above (eq. [29]) is obtained in the Appendix by setting the right-hand side of 
equation (46) zero. 
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Within an order-unity factor, k 2 e^ ~ (w\) 2 when A = 1/k. Recall from §3 that w\ is the 
"typical" value of on lengthscale A; more precisely, it can be defined as the square root of the 
second-order structure function of . Therefore the steady-state scaling wl oc \( 1+a )/ 2 (eq. [22]) 
corresponds to oc fc~( 3+a ); similarly, oc A^ 1- ")/ 2 corresponds to oc k~^~ a \ 

5.1. Energy and Flux in the Kinetic Equation 

The energy per unit mass in up-waves is (|to^| 2 ) /2 = (l/2)(2-7r) -4 f e 1[ d 2 k. So e 1 is propor- 
tional to the energy per unit d 2 k. Conservation of up- wave energy implies that (d/dt) J (ke^)dk = 0, 
assuming isotropy and neglecting dissipation. This can be seen immediately from equation (46) if 
we re- write it as follows (without the diffusive term), 

/•oo 

{d/dt)(ke\) = (A/v A ) / (4 2 - e{)Sl(k,k 2 )dk 2 , (48) 
J o 

since is symmetric, 

/>2tt 

S l (k,k 2 ) = {kk 2 f / sin 2 9 cos 2 9 e[ fcf 2 = S l (k 2 ,k) , (49) 
Jo 1 

and el — e\ is antisymmetric. In comparing our results with those of Galtier et al. (2000), note 
that these authors use the energy per unit k, which they denote E + , i.e., E + ~ ke\ within a 
multiplicative constant. 

The energy flux is the net rate at which energy flows across a given wavenumber k. We 
define it as the integral of the right-hand side of equation (46) over d 2 fe/(2vr), from some particular 
k to infinity, so that the kinetic equation can be written as 

(d/dt){ke\) = -{d/dk)e^ . (50) 

Explicitly, we can write the flux in the form 

£ T= / dp dq(el-el)S l (p,q) . (51) 

JO Jk 

It is trivial to verify that d/dk of this expression is equal to minus the right-hand side of equation 
(46). The positive term in equation (51) is the rate at which energy flows from wavenumbers smaller 
than k to those greater than k; the negative term is similar, but with the origin and destination of 
the energy switched. 

6. Numerical Simulations of Kinetic Equations 

In this section we present the results of numerical simulations which verify our previous heuris- 
tic discussions: in particular, the pinning of spectra and the scaling of the spectra in steady state. 
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It is much faster to simulate kinetic equations than the full equations of motion for and w^. 
There are a number of reasons for this. First, the kinetic equations are only one dimensional, since 
homogeneity has been assumed parallel to the mean magnetic field, and isotropy has been assumed 
in the plane transverse to the mean magnetic field. Second, and more importantly, the averaged 
energies and e 1 are much smoother functions of k than and w^. Thus, a logarithmically- 
spaced grid can be used, which greatly reduces the number of variables that need to be evolved. 
The result is an enormous reduction in computational time. A typical kinetic simulation takes a 
few hours on a PC to reach steady state. A comparable fully three-dimensional MHD simulation, 
would require many months, if not years, on the fastest supercomputers. 

Galtier et al. (2000) perform numerical simulations of the kinetic equation. However, their 
investigation of the imbalanced cascade is incomplete. In particular, they only plot spectra of the 
product e^eK They do not discuss the pinning of the spectra at the dissipation scale, which is 
crucial to the evolution of the cascade. 

For our numerical simulations of the kinetic equations (eq. [46] and the analogous equation) , 
we set 

A * (52) 

This corresponds to absorbing 2A/va into and e^. Our method of integration is very similar to 
that of Galtier et al. (2000). We change variables in equation (46) to give 



8 , r °° 



dt 



k 



r°° -i r\k+k 2 \ pUki t) 

t) = k dk 2 k\ e T (ife2,t)-e T (M) / dh sin 9 cos 2 9 \ 1; ; - vk 2 e\k, t) , (53) 
Jo - J J\k-k 2 \ fe l 



where 9 is a function of k\ and k2 (given in eq. [47]). We evaluate all functions of k on a fixed, 
logarithmically-spaced grid with k = 2*/ 8 , i = 1, 100. At the outer scale, 

k out = 2 1 / 8 = 1.09 , (54) 

and the maximum k is 

AW = 2 100/8 = 57 93 . (55) 

The double integral over k\ and &2 is performed by summing the values of the integrand evaluated 
on the /c-space grid. The factor sin 9 cos 2 9 is averaged in the vicinity of each grid point, i.e., at 
each (/ci,/c2)- Since 9 is a function oik\/k and A^/fc, and since the grid is logarithmic, the averaged 
angular factor can be precomputed and stored in a two-dimensional matrix, each element of which is 
the average in the vicinity of (k\/k, A2/A). 13 We integrate in time with second-order Runge-Kutta. 



13 In evaluating the integral near k out , the values of e T (k 2 < k out ) and e^(fci < fc ou t) are required. Since these 
values are off of the grid, some method of extrapolation is needed. In the runs that we present in this paper, we 
extrapolate with e T , e 1 oc k~ 2 . With this extrapolation, there is no energy transfer from modes with k < k out to those 
with k > fcout; steeper extrapolations would transfer energy. We have also experimented with a natter spectrum for 
k < k out : e T = e T (fc out ), = e^fcout). While this changes the behaviour near k ~ fc out — in particular, it leads to a 
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For the discussions of the simulations that follow, recall that the energies of the waves at a 
given lengthscale A = 1/k are, within unimportant multiplicative constants, 

(w{) 2 ~ fcV(fc) , (^) 2 ~ k 2 e l (k) . (56) 

The cascade times of the up- and down-going waves are, respectively, 

tL(k) ~ (fcV(fc))- 1 , t-Uk) ~ (fcV(fc))" 1 , (57) 

(see eq. [16]). The dissipation time is 

tdUk) = (vk 2 )- 1 , (58) 
(see eq. [24]). In the simulations in this section, the viscosity is 

u = 3 • 10" 5 ; (59) 

this implies that tdiss(femax) = 0.001. 

6.1. Fixed Energies at the Outer Scale 

For our first simulation, we fix e^(k ou t) = 1 and e-'-(fcout) =0.1 throughout the simulation. We 
do this by adding an injection term at k = k out to the right-hand side of the kinetic equation (eq. 
[53]). The two injection terms, ej n - and e^, are adjusted to keep e 1 * (k out ) and e^(/c out ) fixed. At 
the initial time, we set oc k~ 3 and oc k~ 3 (Fig. la). These initial spectra are seemingly valid 
solutions of the steady-state flux relations (eqs. [29] and [30]), with a = and = But the 
spectra are not pinned. When they are evolved in time, it seen that they become pinned to each 
other at the dissipation scale. This pinning happens very quickly — at the dissipation timescale (Fig 
lb). The reason for this pinning was discussed in §3.3, and can be traced in Figures la-c. 

The entire e 1 spectrum adjusts to e} on the timescale tcas(fe ou t) ~ 1- The e} spectrum takes 
longer to adjust, since t<L,s(fc ut) ~ 10- By t = 50, steady state is reached (Fig Id). We discuss 
the resulting steady state spectra in §6.3. We find that, in steady state, the energy injection rates 
required to keep the outer scale energies fixed are 

4 Q . = 0.136 , 4 - = 0.114 ; (60) 

these values are used in our second simulation. 



sharp drop in e T and e 1 between k = fc out and k = 2 1 / 8 fc out by a factor of around 3 — the remainder of the spectrum 
for k > 2 1 / 8 fc out is nearly unaffected. With the k~ 2 extrapolation there is also a drop at k be seen in 

the top panel of Figure Id, for example; but this drop is much less drastic than with the flat spectrum extrapolation. 
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Fig. 1.- 



Simulation of Kinetic Equations with Fixed Energies at the Outer Scale 



-20- 



6.2. Fixed Fluxes 

In our second simulation, we inject energy at fixed rates at k = k out , allowing e^(/c ou t) and 
eH^out) free to evolve. We use the same injection rates as we found in steady state in the previous 
simulation (eq. [60]). We add to (k out ) and e^(k out ) at these fixed rates throughout the present 
simulation. Since we add to e^(/c ou t) and e^(/c ou t) at the same rate as we did in steady state in the 
previous simulation, we expect that our second simulation will reach the same steady state as did 
the first one. We shall show that it does. Note that, strictly speaking, ej n - and are not precisely 
equal to the fluxes. Nonetheless, we ignore this subtlety and refer to the present simulation as one 
of constant fluxes. We shall evaluate the fluxes more precisely below. 

For our initial condition, we use the spectra obtained in steady state in the first simulation, but 
multiplied by a constant; specifically, e^(k) — > e^(k)/10 and e^(k) — > e^(k) ■ 10 (see Fig. 2a). With 
these initial spectra, and with the aforementioned fluxes, the flux relations (eqs. [29] and [30]) are 
satisfied — since they were satisfied before the multiplication and division by 10. But, once again, 
this is not a valid steady state solution because the spectra are not pinned at the dissipation scale. 
When the spectra are evolved in time, the spectra are first pinned (Figs. 2a-c). In this simulation, 
the spectra are initially pinned where teas ~ t^ ss , which is at k ~ 1,000. Therefore the pinning 
timescale is ~ tdiss(^ = 1 ; 000) ~ 0.03. 

As expected, the same steady state is reached as in the first simulation (Fig. 2d). In reaching 
steady state, the two spectra must cross. This happens at t ~ 40, when the two spectra are nearly 
identical, with slopes equal to —3, and with outer scale energy equal to ~ 0.3. The time to reach 
steady state is considerably longer in the constant flux simulation than in the constant energy one. 
Steady state is reached at t ~ 250, as opposed to t ~ 50 in the previous simulation. 



6.3. Energies and Fluxes in Steady State 

We can quantitatively compare the behaviour in steady state with our calculations in §4. 
We consider the steady state reached in the two simulations discussed above. Recall that both 
simulations reach the same steady state, at which point the spectra are as shown in Figures Id and 
2d. For our comparisons, we need the following quantities: v = 3 ■ 10~ 5 (eq. [59]), k out = 1.09 (eq. 
[54]), A/v A = 1/2 (eq. [52]), e^{k out ) = 0.63, and e^k out ) = 0.066. 14 

From equation (27), the dissipation wavenumber is 

fediss = — (feout 6 e T (fe out )eHfc O ut)) 1/2 = 4,400 , (61) 



14 For our first simulation, we actually fixed e^(fc ou t) = 1 and e^fcout) = 0.1. However, as we discuss in footnote 13, 
there is a drop in energy between the first and second grid point. The values c T (fc out ) = 0.63, and e^fcout) = 0.066 
are found by taking the power laws seen in steady state for e T (fc > fc out ) and e^(k > fe ou t), and extrapolating them 

to k — kout • 
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Fig. 2.- 



Simulation of Kinetic Equations with Fixed Fluxes 
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in agreement with the value seen in Figure Id. 
Equation (28) gives 



ln[e T (fc out )/e' l (fco U t)] 
ln[/c out 4 eT (£; ou t)e-L (k out ) (A/ vv a) 2 ] 



a = — TT71 \~7T~/ \oT = °- 14 . ( 62 ) 



which implies that the steady state slopes should be 

din e 1 * 



dink 

din 

dink 



3 + a = 3.14 (63) 
3 - a = 2.86 . (64) 



This is in agreement with the top panel of Figure Id, which shows that the slope of in the 
inertial range, while varying slightly with k, mostly remains between 3.1 and 3.15; the slope of 
is between 2.85 and 2.9. 

To derive the steady-state formulae that we have used thus far in this section, heuristic argu- 
ments suffice. However, to relate energy spectra to fluxes, the kinetic equation is required. From 
equation (33), we should have 

eV = [1.87 • k out 6 e^(k out )e l (k out )(A/v A )} 2 = 0.0043 , (65) 

where the numerical factor 1.87 is obtained from the kinetic equation. A more important application 
of the kinetic equation is to derive the ratio of the fluxes, for which heuristic arguments are useless. 
From equation (32), we should have 

- 1 = a = 0.14 , (66) 

when we use the value of a predicted in equation (62). To compare the above predictions for the 
fluxes (eqs. [65] and [66]) with the fluxes seen in the numerical simulation, we plot the latter in 
Figure 3. We calculated with equation (51), using the steady state spectra shown in Figure Id. 
The quantity S^(p, q) that appears in equation (51) is given in equation (49). Since the kinetic code 
calculates S^(p, q) in order to evolve the kinetic equations (see eq. [53]), it is trivial to modify the 
code so that it can be used to evaluate the flux. To calculate e^, we used the analogue of equation 
(51), with | and j switched. We see in Figure 3 that the fluxes are nearly independent of k in the 
inertial range, with ~ 0.070 and ~ 0.061. These values give e^/e 1 - 1 ~ 0.15, in agreement 
with equation (66). They also give ~ 0.0043, in agreement with equation (65). 

From Figure 3, dissipation affects and for k > 100. Thus dissipation has a big reach. If 
we add the dissipation term to the flux-conservation form of the kinetic equation (eq. [50]), then 
we see that in steady state the kinetic equation must satisfy 



(67) 
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Fig. 3. — Steady State Fluxes and Flux Gradients 

Fluxes and flux gradients are shown as bold solid and dotted lines. Also shown as unlabelled thin 
lines that nearly overlap the flux gradients are uk 4 e^ (thin dotted line) and vk A e^ (thin solid line). 
The spectra from which the fluxes were extracted are shown in Figure Id. 
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The left-hand side of this equation gives the rate of energy increase (per logarithmic band in k) 
due to the /c-space divergence of the flux; in steady state, it must be balanced by the energy loss 
rate due to dissipation. In Figure 3, we plot — (d/d\n k)e^ as a bold dotted line; we evaluated it 
by taking the numerical derivative of the displayed e) . The quantity uk 4 e^ is the unlabelled thin 
dotted line that nearly overlaps — (d/dln k)e^ . Thus equation (67) is satisfied in the simulation. 
The same is true for the down waves, which are shown as solid lines in Figure 3. 

6.4. Decaying Turbulence 

For our third simulation, we allow the spectra to decay without injecting any energy. The 
initial condition is the steady state spectra from the previous simulation, see Figure 4. At the outer 
scale, the cascade time of the down-waves is ~ l/e^(k ou t) ~ 1, so e^(k ou t) decays on this timescale. 
More precisely, from the values plotted in Figure 4, oc exp(— 0.5t) at fixed k. 

The cascade time of the up-waves is much longer. Initially, iLs ~ 10 at the outer scale; as 
decays, iJas increases exponentially. So (k out ) does not evolve. At larger fc's, e} is cut off by 
dissipation; the dissipation wavenumber is ~ {vt)" 1 ^ 2 ~ 200i~ 1 / 2 . 

The spectra appear to evolve in a self-similar manner, with the spectra remaining pinned at 
their dissipation scale. The end result is that the energy in the down-waves disappears, while 
the energy in the up-waves is nearly unchanged. Decaying weak turbulence is unstable: an initial 
imbalance between up and down waves is magnified exponentially. A similar instability occurs in 
strong turbulence, as suggested by Dobrowolny et al. (1980) in the context of the solar wind, and 
as seen in numerical simulations of strong turbulence (Maron & Goldreich 2001, and Cho et al. 
2002). 

7. A Model of the Kinetic Equations: Coupled Diffusion Equations 

It is instructive to model the kinetic equations with two coupled diffusion equations. Inves- 
tigation of these model equations illustrates that the pinning of spectra is quite general; the only 
requirement is that the cascade time of one type of wave be inversely related to the energy of the 
other type. A second reason for considering model equations is that they can be simulated much 
faster and more easily than the kinetic equations. Finally, by contrasting the kinetic equations with 
the model equations, we can gain insight into the behaviour of the kinetic equations. 

7.1. Derivation of Coupled Diffusion Equations 



Our model equation is 



d | _ d _j 

~dt = ~dk e ' 



(68) 




Fig. 4. — Decaying Simulation of Kinetic Equations 
Spectra of e 1 * (dotted lines) and e 1 (solid lines) at the following times: t = 0, 25, 50, 75, 100, 125 

(from top to bottom). 
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where ke' is the energy per unit k, and is the flux. This equation is similar in form to equation 
(50). Recall from §5.1 that we define the flux as the energy flow rate across a fixed k = \k\. It is the 
"1-D flux," since we average over angles; similarly, ke} is the "1-D energy density." For e\ we choose 
a form that depends on the following quantities evaluated at k: e^, e\ and due 1 ■ By contrast, the 
kinetic equations depend on and evaluated at a range of wavenumbers, so unlike the kinetic 
equations, the diffusion equations are exactly local in k. For the cascade time of the energy in 
the up waves to have the correct form (eq. [16]), i.e., tLs = (X/w^) 2 (va/ 'A) ~ (/c 4 e^) _1 (f a/A), we 
choose 

e ~T = _A e i fc 6-/3!_/y+i e A (69) 

va ok\ J 

where (3 remains to be specified. Of course the above relation for the flux gives the correct steady- 
state scaling: eW oc k~ 6 . The evolution equations for are the same as for (eqs. [68] and [69]), 
with |'s and j's interchanged. 

With e 1 oc k~( 3+a \ oc k~( 3 ~ a \ we can relate the steady-state fluxes to the energies: 



gT = f( a )[k 6 eUk)e l (k)]— , (70) 

va 

& = f(-a)[k 6 e ] (k)e l (k)}— , (71) 

VA 

where 

f(a) = 2 + a - (3 . (72) 

These relations are analogous to equations (29) and (30). To make the analogy closer, we choose (3 
so that / has the same dependence on a as does / in the limit that the cascade is nearly balanced 
(i.e., \a\ <C 1). Since e^/e^ = f(a)/f(—a) ~ 1 + a/(l — (3/2) in this limit, we see by comparison 
with equation (32) that we should choose (3 = 0. 

Collecting results, our model equations are two coupled diffusion equations: 



A 



e 



A* * (V]fe 6 ^]feeT)-i/fcV, (73) 



d t 

— e . , , 

dt VAkdk\ dk 



-e* 



dt va k dk 

after including dissipative terms. 



(e^k 6 — ke l ) - z/fcV , (74) 
V ok J 



These equations, while much simpler than the kinetic equations, share many of the same 
features. Since the cascade time of the up-going waves is inversely proportional to the energy of the 
down-going waves (and vice versa) , the argument for the pinning of the spectra in weak turbulence 
(see §3.3) applies here as well. Moreover, in steady state, these diffusion equations suffer from the 
same degeneracy as does the kinetic equation: the constancy of e 1 ' and is insufficient to determine 
the scaling of e 1 ' or separately. This degeneracy is partially broken by the dependence of the 
fluxes on the slopes (eqs. [70] and [71]); we have chosen this dependence so that it is the same as 
for the kinetic equations in the limit that the cascade is nearly balanced (\a\ <C 1). In the following 
section, we present numerical simulations of these equations. 
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7.2. Numerical Simulations of Coupled Diffusion Equations 

We run two simulations of the coupled diffusion equations. As with the kinetic simulations 
described in §6, the first simulation has fixed energy at the outer scale, and the second simulation 
has fixed flux. And, as before, 2k/ v a is set to unity, functions of k are evaluated at k = 2*/ 8 , 
i = 1, 100, and the viscosity is v = 3 ■ 10~ 5 . 

7.2.1. Fixed Energies at the Outer Scale 

For our first simulation, we fix e^(fc ou t) = 0.63 and e^(k ou t) = 0.066. 15 The evolution is shown 
in Figure 5. It is very similar to the evolution of the kinetic simulation (Fig. 1). As before, 
the spectra are pinned quickly. And since, by design, the predicted steady state relations for the 
diffusion equations are the same as those for the kinetic equations (eqs. [61]-[66]), the steady 
state behaviour of the two simulations are nearly identical; compare Figures 5d and Id. There are 
two differences worthy of note. First, the slopes of e} and do not have a spike near k = fc out . 
The presence of such a spike in the kinetic simulation is due to the extrapolation of the spectra to 
k < fcout (see footnote 13). Since the diffusion equations are exactly local in k, such an extrapolation 
is unnecessary, and so the behaviour is much smoother near k = k ut- Second, in steady state, the 
diffusion simulations yield spectral slopes that change more gradually as a function of k near the 
dissipation scale; compare the top panel of Figure Id with that of Figure 5d. 

7.2.2. Fixed Fluxes 

For the second simulation, we inject energy at k = k out at a fixed rate, equal to that found in 
steady state in the previous simulation. For the initial condition, we divide the previously found 
steady state value of by 10, and we multiply e 1 by 10. The evolution is shown in Figure 6. It is 
very similar to that found in the corresponding kinetic simulation (Figure 2). 

8. Large Magnetic Prandtl Number 

In incompressible MHD, diffusive dissipation of kinetic and magnetic energies is accounted 
for by addition of the terms vV 2 v and 7]V 2 B to the right-hand sides of equations (1) and (2), 
respectively, where v is the viscosity and n is the resistivity. The magnetic Prandtl number is 
defined as 

Pr = u/rj . (75) 



J We choose these numbers so that the spectra will have the same amplitude here as in the kinetic simulation; see 
footnote 14. 
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Fig. 5.- 



Simulation of Diffusion Equations with Fixed Energies at the Outer Scale 
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Fig. 6.— 



Simulation of Diffusion Equations with Fixed Fluxes 
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Until now, we have assumed that Pr= 1. In this section, we examine weak MHD turbulence when 
Pr> l. 16 

With arbitrary v and rj, the equation for (eq. [4]) becomes, 

d t w^ +vaO z w^ = -w l -VwT -VP + ^(i/ + tj)VV + -{v - r])X7 2 w l ; (76) 

the equation for is the same but with f and [ interchanged, and va — ► — va- To derive the 
kinetic equation for e\ we first Fourier transform the above equation; the transformed dissipation 
terms are — (k 2 /2)(u + vi)w\ — {k 2 /2){v — r])w^ k . When we take the dot product of the Fourier- 
transformed equation with wl, and then take the expectation value of the result, we can neglect 
the dissipation term proportional to (w^-wj^, since it is much smaller than the term proportional 

to (^(iVfc) 2 ^. This is because in upgoing waves at fixed z^ = z — v^t are uncorrelated with the 
downgoing ones (see §A.l). As a result, the kinetic equation for high magnetic Prandtl number is 
the same as that for Pr= 1 (eq. [46]), but with uk 2 e^ replaced by {y + r])k 2 e} /2 ~ vk 2 e} /2. This 
might appear to be a surprising result. In the limit of vanishing rj, the energy spectra are cut off 
below the viscous scale, i.e., the scale set by v (eq. [27], within a factor of 2). One might have 
expected that the viscous scale should only affect the kinetic energy spectrum. Yet the magnetic 
energy spectrum is also cut off below the viscous scale. This can be understood as follows: upgoing 
waves each have equal magnetic and kinetic energies, and similarly for downgoing waves. As an 
up-wave with lengthscale slightly larger than the viscous scale is gradually cascaded to sub-viscous 
scales, its kinetic energy is dissipated by viscosity. Since the cascade time in weak turbulence is 
much longer than the waveperiod, as the up-wave's kinetic energy is dissipated, its magnetic energy 
is converted into kinetic energy so that its magnetic and kinetic energies can remain nearly equal. 
As a result, both kinetic and magnetic energies are dissipated at the viscous scale. 



9. Hyperviscosity and the Bottleneck Effect 

In many different types of turbulent cascades, the energy spectrum exhibits a hump on scales 
slightly larger than the dissipation scale (e.g., Falkovich 1994). The hump is particularly pronounced 
in numerical simulations that use "hyperviscosity," a trick whereby the diffusive term vk 2 v is 
replaced by u n k n v, where n is typically 4 or 8, and v n is the hyperviscosity; the resistivity is 
modified in a similar manner (e.g., Borue & Orszag 1995, Biskamp 2000). With this trick, a 
smaller part of the spectrum is affected by dissipation, so a longer inertial range can be simulated 
with a fixed resolution. Although this trick does work, there is a problem: the spectrum on scales 
slightly larger than the dissipation scales is made flatter. The energy, in effect, is backed up. This is 
the bottleneck effect. It can be particularly problematic in simulations of strong MHD turbulence, 



Strong MHD turbulence with PrS> 1 occurs in many astrophysical settings. 
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where the energy backup can affect lengthscales considerably larger than the dissipative scales 
(Biskamp k Miiller 2000). 

This motivates us to investigate the bottleneck effect in weak turbulence. If the bottleneck 
effect appears, its interpretation will be simpler than in strong turbulence. We perform two simu- 
lations of the kinetic equation, one with hyperviscosity of the form v^k A , and the other with vgk 8 . 
In each simulation, we fix the energies at the outer scale, e^(fe out ) = 1, e^(A; out ) = 0.1, and allow the 
spectra to reach steady state. The steady state spectra of /c 6 eW are shown in Figure 7, offset for 
clarity. Also shown is the simulation with ordinary viscosity described in §6.1. From this figure it 
is apparent that weak turbulence suffers from the bottleneck effect; the effect becomes larger with 
increasing hyperviscous exponent. 

The bottleneck effect can be understood as follows. Consider an up-wave with a lengthscale 
slightly larger than the dissipation scale. It is cascaded by down-waves that have slightly different 
lengthscales than its own. Hyperviscosity gives a sharper dissipative cutoff to the down-wave 
spectrum than ordinary viscosity. Therefore a hyperviscous simulation has fewer down-waves in 
the vicinity of the dissipation scale, and the cascade time of the up-wave is longer. For the up-wave 
energy flux to be independent of lengthscale, a longer cascade time implies a larger energy. As a 
result, the spectrum is flatter on scales slightly larger than dissipative ones. Falkovich (1994) offers 
a similar explanation for the bottleneck effect in hydro turbulence. 

10. Discussion 

In this paper we describe imbalanced weak turbulence and solve the steady state cascade. In 
a future paper, we will extend the result to the strong cascade. One of our ultimate goals is to 
develop a theory of imbalanced strong turbulence to apply to the solar wind, where imbalance is 
observed. 

Although strong turbulence is more generally applicable than weak, the latter is a simple and 
illuminating model. There are a number of issues in strong turbulence that are not understood. 
Weak turbulence can be used as the first step in explaining them. For example, in this paper we 
examined the effect of a large magnetic Prandtl number on weak MHD turbulence. We also found 
that the bottleneck effect appears in weak turbulence, where its interpretation is straightforward. 
There are a number of other issues that we intend to investigate in weak turbulence as a prelude 
to understanding them in strong turbulence; for example, reconnection and the turbulent dynamo. 

Research reported in this paper was supported by NSF grant AST-0098301. 
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Fig. 7.— The Bottleneck Effect 

Three spectra, offset for clarity; the top spectrum has viscous and resistive terms vsk s , the middle 

has V4 k 4 , and the bottom has vk 2 . 
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A. Appendix: Kinetic Equation in Weak Turbulence 
A.l. Preliminaries 

Kinetic equations describing weak MHD turbulence are derived by Galtier et al. (2000, 2002). 
In this appendix, we present a more physical derivation. We compare the two derivations in §A.5. 

We consider the evolution of in a plane that is transverse to z and moving with velocity 
vaz, i.e., with fixed = z — v^t- Recall that £ is in the direction of the mean magnetic field. 
Changing variables from z to z^ in equation (8) gives 



d_ 

dt 



x,y,z< 



= -w i -V±w^ + V±V ± 2 ('V±w i :V±w^) . (Al) 



In weak turbulence, the parallel length over which disturbances are correlated, A, is small; in 
particular, in the time that an up- wave crosses a down-going wavepacket, At ~ A/va, its distortion 
is less than unity: Aw\/w\ ~ (w^/va)A/X <C 1 (eq. [15]). 17 Therefore, undergoes small 
uncorrelated changes each At ~ A/va, and its evolution is analogous to a random walk, with step- 
size ~ w^(w^/va)A/\ <C 1. Our goal in this appendix is to quantify the random walk; the 
resulting equation is the kinetic equation. The kinetic equation is only valid when the turbulence 
is weak: its derivation hinges on the assumption that the time for up-waves to cascade is longer 
than the correlation time of down-waves at fixed z\ i.e., that (w^/va)A/\ <C 1. 

Equation (Al) is linear in . Nonlinearity arises because modifies (through eq. [10]); 
this modification backreacts on through equation (Al). Nonetheless, this backreaction can be 
neglected when deriving the kinetic equation. The reason is as follows. Every At ~ A/va, the 
backreaction changes w\ by ~ w^ x [(w^/vA)A/X][(wj i /vA)A/X], which is smaller than Aw| by the 
small factor (w\/va)^/\ ^ 1. Furthermore, since this backreaction is uncorrelated on timescales 
larger than A/va, it only effects a small change in the step-size of the random walk. 

Therefore, in deriving the kinetic equation, we interpret equation (Al) as a linear equation 
for w;1 (at fixed z 1 * ); at fixed z 1 * is viewed as a function with known statistical properties. To 
simplify our derivation, we first solve a model problem: the linear random oscillator. The extension 
to weak turbulence will then be straightforward. 



A. 2. A Toy Problem: the Linear Random Oscillator 

We consider the evolution of a simple random oscillator ip: 

j t m = iA(t)m , (A2) 



7 We neglect the factor of 2 associated with the fact that the relative velocity of up and down waves is 2va- 
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where ip is a complex scalar and A is a real random variable; the factor i ensures that the energy 
| | 2 is conserved. 18 Our goal is to calculate the evolution of (ip(t)), where angled brackets denote 
an average over an ensemble of A's, not over time. We assume that the values of A at two different 
times are statistically independent of each other whenever the two times are separated by more 
than the correlation time, r corr ; we make this more precise in footnote 20. For simplicity, we take 
A to have zero mean, (A(t)) = 0; the extension to A with non-zero mean is trivial. We assume that 
the statistical properties of A — such as T cor r and A ims = ( J 4 2 ) 1 ^ 2 — vary only on timescales much 
larger than r CO rr- Note that ip corresponds to in equation (Al); A corresponds to w^, or, more 
specifically, to w^/X; and r corr corresponds to A/va- 

There are two limiting regimes for the random oscillator, depending on whether Arms r cor r is 
less than or greater than unity. In the following, W6 tclkc -Arms Tcorr "C 1; this regime corresponds 
to weak turbulence. The change in ip within the time T cor r is of order A Tms T COTT ip , which is much 
smaller than ip. Thus the correlation time of A is smaller than the "cascade time" of ip. Since A 
is uncorrelated on timescales larger than T CO rr> ip undergoes small uncorrelated changes each T cor r) 
and thus its long-time evolution is a random walk. Intuitively, we expect that the time evolution of 
the statistical properties of ip — in particular, (ip) — can be represented by a differential equation. 
Our goal in this section is to derive the differential equation, and to understand the approximations 
that are made in deriving it. 

The solution of equation (A2) is simply 

= ip(0) • exp(iA) , (A3) 

where A 1 = J* A(t')dt'. We expand as follows, 

i/>(t) = V(0) • (1 + %A t - (AAj + ...), (A4) 

where (A4*)* = ft A(t') f£ A(t")dt"dt'; note that (1/2)(A) 2 = (A4*)* is an identity for any A(t). 
Equation (A4) may be thought of as an expansion in powers of A; it is only valid for short times 
after t = 0. 

We evaluate (ip(t)) as follows. Since A(t) is unaffected by ip(0), A(t) and ip(0) are uncorre- 
lated. 19 So 

W)) = W0)> • (1 + i (A*) - ((AA*)*) + (. . .)) . (A5) 

Since (A) = 0, we have (A*) = 0. Next, we consider ((AA*)*). The values of A at two different 
times "statistically overlap" only when the two times are separated by less than the correlation 

18 For a textbook discussion of the linear random oscillator, see van Kampen (1992). 

19 In fact, A(t > 0) and ip(0) are slightly correlated: ip(Q) is affected by A(ti < 0), which is in turn correlated with 
A(t2 > 0) as long as t 2 — ti < r corr . Nonetheless, in the limit A rms r corr <C 1 that we are considering, tp only changes 
by a small amount in the time r corr . So the correlation between ^4(0 < t < t coii ) and rp(0) is unimportant for the 
evolution of V'W- 
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time; i.e., (A(t')A(t' + At)} is non-zero only if At < T cor r- Thus (AA 1 ) ~ A 2 ms T corr when t > r corr , 
and 

((^)') = 4wt. a (A6) 

Therefore 

^ (V>(0)) (1 - A 2 rms r COII t) , r corr < t < r corr /(^ rms r corr ) 2 . (A7) 

Equivalently, 

— (-0(f)) = _^2 mgTcorr ^ assuming A rms T cor r < 1 . (A8) 

This equation is the main result of this section. Van Kampen (1992) gives a a more mathematically 
rigorous derivation of it than we do. Equation (A8) can be understood as follows: since ip changes 
by Atp ~ ^4rmsTcorrV' m the time T CO rr> then after N steps, each r cor r long, the change in tp is -\/N Aip. 
Thus for order unity changes in tp, (tp/Aip) 2 steps are required. The resulting time, the "cascade 
time," is 

r cas = (V>/AV0 2 r corr ~ l/(^ msTcorr ) , (A9) 

as in equation (A8), and T COTT /T CSiS ~ (^4 r ms' r corr) 2 "C 1. Although (ip(t)) decays to zero on the 
timescale r cas , the energy \tp\ 2 remains constant; the "cascade" of (tp) is analogous to phase mixing. 

We may now proceed to derive the kinetic equation in weak turbulence. Before doing so, 
we consider the linear random oscillator in more detail. If the reader is satisfied with the above 
derivation of equation (A8), the following subsection may be skipped. 



A. 2.1. The Validity of Perturbation Theory and Goldreich & Sridhar (1997) 

Goldreich & Sridhar (1997) incorrectly claim that perturbation theory fails in weak turbulence, 
i.e., in their "intermediate" turbulence. 21 We explain their claim and its resolution in the context of 
the linear random oscillator. In the process, we clarify the validity of the perturbation expansion. 

We consider the terms neglected in equation (A7). For example, the fourth-order term is 
(ip^) I {ip(0)) = ( K {A[A(AA t ) t ] t } t y In this quadruple time integral, A is evaluated at four different 
times. Whenever two of these times are separated by less than r CO rr; the values of A at these two 
times "statistically overlap," and hence can give a non-zero contribution to the total integral. If 
two of the times are separated by less than r cor r; and the other two times are also separated by 



20 For this equation to be approximately valid, A must decorrelate sufficiently rapidly that {A(t')A(t' + At)) goes 
to zero faster than 1/Ai for large At > r COTT ; otherwise, ((^4^*)') rises faster than the first power of t, and our 
derivation is invalid. 

21 A note on terminology: the turbulence that we and that Galtier et al. (2000) call "weak," Goldreich & Sridhar 
(1997) call "intermediate." They call it "intermediate" precisely because of their claim that perturbation theory is 
invalid. 
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less than T cor r, then the integrand does not vanish, even if the first two times are separated from 
the second two times by more than T CO rr* 

So the quadruple time integral yields approximately 

<^)> / <V(0)> ~ (^ 2 ms W) 2 = (t/r cas f. 

Since (ip^) / (V>(0)) = —A 2 ms T corr t = — t/r cas , the contribution of the fourth-order term is as 
large as the second-order term after the time t = r cas . Similarly, if we consider the contribution to 
^( 2n )^ / (^(O)) from correlating pairs of A, the result is ~ (i/r cas ) n ; so all terms are of comparable 
value when t = r cas , and it seems that the lowest order term is inadequate. This, in effect, is the 
claim that Goldreich & Sridhar (1997) make in the context of weak turbulence. 

It is incorrect; although equation (A7) for (ip(t)) is only valid for t <C r cas , equation (A8) for 
(d/dt) (Y>(i)) is approximately valid for all times, with corrections of order powers of T COTT /T cas <C 
1. Since ip(t) undergoes small uncorrelated changes every timestep of length T cor r, we expect on 
physical grounds that the evolution of its statistical properties should be governed by a differential 
equation that is invariant under time translations. 22 Equation (A8) is the only such equation whose 
small-time behaviour is given by equation (A7). Its right-hand side may be interpreted as the lowest 
term in a perturbative expansion in t coii /t CSlS . All of the contributions to (ip(t)) that are of order 
(t/T cas ) n must be derivable from equation (A8). For example, if ^4 2 ms r cor r = 1/fcas is constant, then 
equation (A8) has the solution (tp(t)) / (ip(0)) = exp(-i/r cas ) = 1 - t/r cas + (l/2)(t/r cas ) 2 + — 
Terms of order (t/r cas ) n in (ip) are generated from the lowest order term. 

We can solve the one-dimensional oscillator exactly when A(t) is Gaussian, and thereby illus- 
trate the validity of equation (A8). If A(t) is a zero- mean Gaussian random variable, then so is A t , 
and (exp(iv4*)} = exp(- ((A 1 ) 2 ) /2). 23 It follows from equation (A3) that 

Wt)> = W0)>-exp(-<(AA')*», (A10) 
assuming that (ip(0)) and A(t) are uncorrelated. Equivalently, 

j t W = - {A A 1 ) <V>> , (All) 

which agrees with equation (A8) when (AA 1 ) = A 2 ms r CO rr- 

Although we are mainly concerned with the limit A rms T corr C 1, we conclude this subsection 
by briefly discussing the opposite limit, A rms T cor r 3> 1- This limit sheds some light on strong 
turbulence, which corresponds to the case ^4rms r corr ~ 1- When A rms r CO rr 3> 1> perturbation theory 
is inapplicable because ip does not undergo small uncorrelated changes every T cor r- Rather, its 
cascade time is ~ l/A vms , which is much shorter than the correlation time of A. Therefore A is 
nearly constant in the time that (if}) cascades. The inapplicability of perturbation theory makes 



It should be invariant on the timescale r CO rr; ^4 rm s and r cor r are effectively constants on this timescale. 

23 Proof: a zero-mean Gaussian x has probability distribution P(x) = (27TXr ms )~ 1 ^ 2 exp(— a; 2 /2x 2 ms ); so (exp(ia;)} = 
I-oo P ( x ) exp(ix)dx = exp(-a; 2 ms /2). 
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strong turbulence difficult, if not impossible, to solve. Nonetheless, we can solve the one-dimensional 
oscillator in the "strong" limit. Equations (A10) and (All) are still valid, but (AA*} ~ ^4 2 ms^> 
so {ip(t)} ~ (^(O)) exp(— yl 2 ms i 2 ), and equation (All) is not invariant under time translations. 
The time t = is special because of our assumption that tp(0) and A(t) are uncorrelated. This 
assumption, while innocuous for the "weak" oscillator, is crucial for the "strong" one. 



A. 3. Derivation of the Kinetic Equation 

We derive the kinetic equation in Fourier space. We Fourier transform equation (Al) in x and 
y (but not z), denoting the transform of by wj^: 

w k = J d2x ^ w ^ exp(-ifc • £Cj_) , (A12) 
where k is purely transverse (k z = 0); similarly, wj^ is the Fourier transform of 

Since is perpendicular to both k and z, it only represents a single degree of freedom. Hence 
it is convenient to define a scalar potential tp^ by wj^ = i(fexl)^|,, with k = k/k. Similarly, 
wj^ = i(kxz)ipfe. The Fourier transform of equation (Al) is then 

d_ 

di 

where 



k Jk^ = j Sp ^fc,p(*)4(*) ' ( Al3 ) 



A k,p(t) = a k,p^ k -pV ( Al4 ) 
-z-{kx p)k ■ p f P . (A15) 



(2vr) 2 v \k-p\ 

We suppress the functional dependences of tp^ and tp^ on z ^ because this equation is evaluated 
at fixed z 1 *; in the following, we replace the partial time derivative by a total derivative, with the 
understanding that z^ is fixed. 24 

We use angled brackets to denote an ensemble average, in a plane with fixed zk We define the 
spectral energy densities e 1 * and e 1 as follows 



™\$W k tt)) = (^ k ml,(t))=eHk,t)5(k + k') (A16) 
w^tyw^t)) = ty i k W l k ,{t)) = e^k,t)5(k + k'), (A17) 



24 Since is real, wj^ = (w where * denotes the complex conjugate, and so ipj^ = j^)*] similarly, 

~ ^ k^*' ^ ^ ouows tnat p ~ ^*-k -p an< ^ a kp = a -k p' ^ e differential energy in up- waves within the fc- 
space area d 2 k is proportional to d 2 k\w^_ | 2 = d 2 fc|?/>^ | 2 . Conservation of up- wave energy necessitates p = —Ap ^ 
and ajf p = ~ a p fc> as can De ver i ne d from equations (A14) and (A15). 
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where and are real; 5(k) is a two-dimensional Dirac delta-function that follows from homo- 
geneity, i.e., from the assumption that (to? (x±) ■ (x± + Ax±)) is independent of x±. 



We shall derive the evolution equation for the bilinear quantity {^jjl^j^'J- F rom equation 
(A13), 

Jt^Uk' ] = J d2 P A k,p[^p\ + ( fe ~ k ') > ( A18 ) 

where (k <-> fc') represents a second term that is the same as the first, but with k and k' inter- 
changed. 

From §A.l, tp^ (and hence ^4^p) can be viewed as evolving independently of tp^ , since the 
backreaction is negligible at fixed z} . Therefore, equation (A18) is nearly identical to that of the 
simple random oscillator (eq. [A2]); [^.V^'] an d ^kp m ^ ne aDOVG equation correspond to tp and 
A, respectively, in equation (A2). Since is a random function with zero mean and correlation 
time Tcorr = A/va, so is Afcp. We solve equation (A18) it in the same manner as we solved the 
random oscillator, using perturbation theory. Expanding in A k p, we find to zeroth order: 

[V^V^'] (0) = constant . (A19) 

To first order: 

[ W>L' ](1) = / d2 p A k,pK'^p ]{0) + (fc ~ k ' ] ' (A20) 

where A^ = J * A k p(t')dt' . To second order: 

^k' ]{2) = J d2 P d2 ^ A k,p A k',q^Wq^ ] + (^4,S)H'^ ](0) } + {k " k ' ] ■ (A21) 

The sum of these last three equations is directly analogous to equation (A4) for the simple oscillator. 
Using the same reasoning here as we did for the oscillator in deriving equation (A8), and assuming 
that the correlation time (A/va) is short, we take the time derivative of the expected value of 
equation (A21); we get, after setting (AA*) = (A/va) (AA) (with the appropriate subscripts on 
A), and after using equation (A16) and integrating out the delta functions on the right-hand side: 

§{k + k')^(k, t) = ±J d 2 p{ (A k ,pA k , t _ p ) e T(p, t) + (A k! pA p _ k ,) J(k>, *)} + (*" *0 ■ 

(A22) 

We re-express this equation by using equations (A14) and (A17) and the relations a k p = a * k p > 
a kp = ~ a pk (^ rom footnote 24), and e^(k,t) = e^(—k,t): 



d 
dt 



e\k,t) = 2— [ d 2 p(e 1 (p,t)-e\k,t))\a k \ 2 e l (k-p,t) . (A23) 
va J 



Since each upgoing plane with fixed interacts with many statistically independent downgoing 
wavepackets before cascading, it is reasonable to assume isotropy in planes transverse to z; after 
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inserting equation (A15) for p and changing variables, 

A f°° r f 27r pUki t) 

e \k,t) = —k 2 dk 2 kl e t (k 2 ,t) - e 1 (k,t)] d9 sin 2 6 cos 2 - , 2 , (A24) 

V A JO ^ JO 



dt 
where 

fci = (A; 2 + k\ - 2kk 2 cos 6) 1/2 . (A25) 

We have redefined A to absorb the factor of Air 3 , i.e., A/(47r 3 ) — > A. Thus far, we have only 
considered a single plane with fixed z 1 * . If we assume that the turbulence is homogeneous in z, the 
z-average of the above equation is trivial since there is no explicit dependence on z. We can simply 
re-interpret angular brackets to denote z-averages in addition to an ensemble average. 

We emphasize that our derivation of the kinetic equation (eq. [A24]) is only valid when 
the correlation time, A/va, is much smaller than the cascade time, (va/ A)(£; 4 ef) _1 . Otherwise, 
equation (A22) does not follow from equation (A21). 



A. 4. Steady State Fluxes 

In steady state, we set the right-hand side of equation (A24) to zero. We substitute power law 
solutions, e^(k) oc fc~( 3+a T), e^(k) oc k~( 3+a ^ , in which case the right-hand side becomes 

eHk )e^k )k^ +ai ^-\k 2 £ dk 2 k 2 3 (k 2 -^) - fc-(3+«T)) (i _ (k 2 /kr +a i) • 

J desm 2 ecos 2 e[k 2 + k 2 2 -2kk 2 cos6) . (A26) 

Note that the square-bracketed term is independent of ko. The above expression follows after 
breaking the k 2 integral into two pieces: one from to k, and the second from k to oo. Then, the 
change of variables k 2 — > k 2 /k 2 is made in the second piece (a "Zakharov transformation"), so that 
its limits are now from to k. Finally, this second piece is combined with the first. 

In steady state, equation (A26) must vanish, so 

«i = -«T • ( A27 ) 

Since w\ ~ (k 2 e^) 1 / 2 oc k~( 1+a i^ 2 , and similarly for w\ , equation (A27) is equivalent to w\w^ oc A 
(eq. [19]). The vanishing of <9 t e^ in steady state yields the same relation as equation (A27), and so 
does not give new information. 

The flux associated with e} is given by integrating the right-hand side of equation (A24) over 
d 2 k/{2ir) from a particular k to infinity. The steady state flux is thus given by integrating equation 
(A26), with ai = — a-f. The result of the flux integration is 

J = f(ai)\e\k )et(k Q )k$—] , (A28) 
l va J 
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where /(a-f) is a dimensionless function of ay. 

1"^ \ — (5— at)/2 

/(a T ) = y dxx 3 lnx(x~ (3+a T) - l)y sin 2 cos 2 0(1 + x 2 - 2x cos 0J , (A29) 

A technical note: although equation (A26) vanishes in steady state, its integral over gives a 
factor of Q| + a i in the denominator, which also vanishes in steady state. This 0/0 ambiguity can 
be resolved by considering the limit as a^ + aj_ — > (L'Hdpital's rule). 

Equation (A28) and the corresponding equation for e^-, are equivalent to those in the heuristic 
discussion (eqs. [29] and [30]). 

Galtier et al. (2000) plot the function f(a^) after numerically integrating the steady-state flux 
integral. 25 They show that — 1 < < 1, that the steeper spectrum always carries more flux (i.e., 
f( a l)/f(~ a l) is a monotonically increasing function of «j), and that in the limit that a^ — > 1, 
e t/ e i _» qq # (Clearly, this also implies that in the limit af — > —1, e^/e^ — ► 0.) In §4.3.1, we explain 
the physical reason why — 1 < a-f < 1. When this condition is violated, the cascade becomes 
nonlocal. Galtier et al. (2000) find infinite fluxes when these inequalities are saturated because 
they consider an infinitely extended spectrum, which leads to unphysical results when the cascade 
is nonlocal. 

As discussed in §4.2, in steady state we are primarily concerned with the case that the up- and 
down-going fluxes are comparable, so \a^\ <C 1. In this limit, we linearize / about a^ = 0, yielding 
approximately 

/(a T ) ~ /(0) • (1 + 0.5a T ) , |a T | « 1 , (A30) 

after numerical integrations of equation (A29). This result is used in §4.2 to calculate the energy 
spectra given the fluxes and Although the value /(0) is much less important than the 
dependence of / on ay we use it when discussing the results of our numerical simulations; it is 
given by 

/(0) ~ 1.87 . (A31) 

A. 5. The Zero-Frequency Mode 

Weak MHD turbulence is often described as being based on three-wave resonances, in which 
one of the three waves has zero frequency (e.g., Shebalin, Matthaeus, & Montgomery 1983, Galtier 
et al. 2000, 2002). However, the interpretation of a zero-frequency wave is unclear. After all, 
it should require an infinite amount of time for three interacting waves to "realize" that one of 
them has zero frequency. In this subsection, we clarify the role of the zero-frequency mode. In the 
process, we shall compare our derivation of the kinetic equation with that of Galtier et al. (2002). 



More precisely, their figure 2 is proportional to [f(oif)f(— a-f)] 1 ^ 4 , and their figure 3 shows /(a-f )//(— a-f ). 
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Note that in our derivation in §A.3 we avoided discussing Fourier wavemodes, since we did not 
Fourier-transform in t or in z. 

The evolution of the up- waves at fixed z 1 * = z — v £t is given in equation (A18). However, for 
clarity, in this subsection we shall consider the random oscillator instead (eq. [A2]): 

dip/dt = iAvp . (A32) 

It is a simple matter to extend our discussion to weak MHD turbulence by replacing tp with V^V^'' 
replacing iA with A^ p , and summing over the appropriate indices. 

In §A.2 we derived the equation for (ip), assuming that the cascade time of (ip) is longer than 
the correlation time of A: 

d(ip) /dt = -(AA t )(iP) , (A33) 
where ^AA*) = J*/_ T (A(t)A(t')) dt' ~ A^ ms T COTT ; the value of T is unimportant as long as T > t coit . 

We can see how the zero- frequency mode enters by re- writing ^AA*) in terms of A w , the Fourier 
transform of A: 

^ - L £ II £ <^> L*"* ■ im 

To compare this result with that of Galtier et al. (2002), we shall use similar notation. We denote 
the power spectrum of A^ by q u : 

AyA.A = q u 5{J + u) ; (A35) 



the Dirac-delta function results from the time-invariance of the statistical properties of A. We also 
define 

A(w) = — e™\\ - e iwT ) . (A36) 
iuj 

With these two definitions, equation (A34) is 

/°° dt i r i 
^^[e-^AH] . (A37) 

If this equation is inserted into equation (A33), the resulting equation is equivalent to equation (7) 
in Galtier et al. (2002). As these authors argue, when T — > oo, exp(— iujt)A(u) — > ttS{u}), so 

/OO J 
^? w <5(a;) , (A38) 

which corresponds to equation (8) in Galtier et al. (2002). From this expression, we can see why 
the zero-frequency mode enters: it is simply a consequence of the term (AA*). It does not require 
an infinite time for (ip) to interact with 9^=0; rather, the only quantity that enters into qoj=o is the 
value of A within the time T cor r of t. We can see this more clearly by explicitly writing the expression 
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for q w . Of course, the power spectrum q u is simply the Fourier-transform of the correlation function 
of A: 



One of our reasons for discussing the Fourier-space picture in detail is that there have been 
a large number of confusing remarks about it in the literature. For example, Galtier et al. (2000, 
2002) claim that, since the kinetic equation is apparently not applicable to the zero- frequency 
mode, there might be a "condensation" of zero-frequency modes. Since the zero-frequency mode 
is so important, this condensation might have dramatic implications for the cascade. However, 
as long as correlation times are shorter than cascade times, the kinetic equation gives a complete 
description of the turbulence, and is applicable to the zero-frequency mode as well, which is simply 
given by g w= o = 4-7T (AA*). For this not to be true, large-time correlations would have to build up 
in A (see eq. [A39]). Since A corresponds to tp^ m weak MHD turbulence, large-time correlations at 
fixed correspond to large-distance correlations along the z-direction. But it is impossible for two 
downgoing wavepackets, initially uncorrelated, to become correlated. This is because the equation 
for their evolution is linear: recall that in weak turbulence the backreaction term is negligible. 
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